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Abstract. We consider a nonlinear damped hyperbolic equation in R"^, 1 < n < 4, 
depending on a positive parameter e. If we set e = 0, this equation reduces to the 
well-known Kolmogorov-Petrovski-Piskunov equation. We remark that, after a change 
of variables, this hyperbolic equation has the same family of one-dimensional travelling 
waves as the KPP equation. Using various energy functionals, we show that, if e > 0, 
these fronts are locally stable under perturbations in appropriate weighted Sobolev 
spaces. Moreover, the decay rate in time of the pertiirbed solutions towards the front of 
minimal speed c = 2 is shown to be polynomial. In the one-dimensional case, if e < 1/4, 
we can apply a Maximum Principle for hyperbolic equations and prove a global stability 
result. We also prove that the decay rate of the perturbated solutions towards the fronts 
is polynomial, for all c > 2. 
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1. Introduction 

We consider the damped hyperbolic equation 

euu{^,t) + ut{^,t) = A^u{C,t) + f{u{^,t)) , (1.1) 

where ^ = (^i, ^2, • • • , ^n) ^ R"^, t G R, e > 0, and / : R — R is a nonlinear map. In the 
one-dimensional case n = 1, equations of the form (1.1) arise as mathematical models 
describing various natural phenomena, like the propagation of voltage along a nonlinear 
transmission line, or the random motion of one-celled organisms [DO]. Here we consider 
the multidimensional case also, and we are interested in situations where the parabolic 
equation obtained by taking the limit e — in (1.1) has a continuous family of travelling 
waves (or fronts) propagating into the unstable state u = 0. Sufficient conditions on the 
nonlinearity / for such a situation to occur are discussed for example in Aronson and 
Weinberger [AW]. For convenience, we restrict ourselves to the typical example of the 
Kolmogorov-Petrovsky-Piskunov equation (KPP), which corresponds to f{u) — u — u'^, 
but more general nonlinearities can be treated by the same methods. 

Existence of travelling wave solutions to damped hyperbolic equations has been 
proved by Hadeler [Ha] in a general context. In our case, this is simply done by setting 

u{^,t) = g{Vl + ec^^i-ct) , (1.2) 

and inserting into (1.1). One obtains for g the differential equation 

g"{x) + cg'{x) + g{x) - g{xf = . (1.3) 

It is well-known [AW] that, for all c > 2, this equation has a front-like solution g{x) 
satisfying g'{x) < for all x G R, liuirc^-oo 9{x) = 1, lim^^^+oo fi'(2;) = 0, and g{x) is 
unique up to a translation in the variable x. Therefore, for all e > 0, Eq.(l.l) has a 
continuous family of travelling waves of the form (1.2) indexed by the parameter c > 2. 
It should be noted that the speed of such a wave is no longer c, but c/\/T+~ec^, a 
quantity which is bounded by 1/v^ as c — > oo. This is of course related to the finite 
propagation speed property of equation (1.1). 

The stability of travelling waves for KPP and similar nonlinear parabolic equations 
has been intensively studied over many years. Early results have been obtained using 
comparison theorems based on the Maximum Principle, see [KPP] , [Fi] , [AW] . Combined 
with probabilistic techniques, these methods give a very detailed description of the basin 
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of attraction of the wave [Bn]. In parallel, a local stability analysis of the front in 
suitable weighted spaces has been initiated by Sattinger [Sa] and continued recently by 
Kirchgassner [Ki], Kapitula [Ka], Bricmont and Kupiainen [BK], Gallay [Ga], Eckmann 
and Wayne [EW], using functional-analytic techniques, renormalization group methods, 
or energy functionals. In particular, the decay rate in time of the perturbations in the 
critical case c = 2 has been investigated [Ki], [BK], [Ga]. Similar results have also been 
obtained for higher dimensional equations [MJ] , and for systems of parabolic equations 
[KR], [RK]. 

The aim of this paper is to extend part of the stability results available for the 
KPP equation to the hyperbolic equation (1.1). In particular, using energy functionals, 
we shall show that the travelling waves (1.2) are locally stable in appropriate function 
spaces for all c > 2 and all e > 0. Moreover, using the Maximum Principle for hyperbolic 
equations, we shall prove a global stability result in the one-dimensional case, provided 
e is sufficiently small. Finally, a decay rate as t — > +00 of the perturbations will be 
obtained if c = 2, or if n = 1 and e < 1/4. 

We now proceed to state our results in a more precise way. Given e > 0, c > 2, we 
go to a moving frame using the change of variables 

u{^, t) = v{Vl + ec^ii - ct, 6, • • • , Cn, t) = v{x, y, t) , (1.4) 

where x = vT+~ec^^i — ct and, if n > 1, y = (^2, • • • , ^n)- The equation for v is 

+vt- 2ecvxt = Vxx + \v + cvx + V - , (1.5) 

and by construction v{x, y, t) = g{x) is a stationary solution of (1.5). As in the parabolic 
case, this solution can only be stable if we restrict ourselves to perturbations which decay 
to zero sufficiently fast as a; — > +00. To achieve this decay, we look for solutions of the 
form v{x, y, t) = g{x) + a{x)w{x, y, t), where a{x) = e~'^^ for some 7 > which will be 
fixed later. Then w{x, y, t) satisfies the equation 

ewtt + {l + 2ecj)wt-2ecwxt = Wxx + \w + {c-2-f)wx + {l-c-f+-f'^-2g)w-aw'^. (1.6) 

Since Eq.(1.6) is of second order in time, we shall rewrite it in the usual way as a first 
order system for the pair {w,Wt), and study the stability of the origin {w,Wt) = (0,0) 
for this system in a space which we now describe. 

Function spaces. For all j G N, we denote by — if-^(R") the usual Sobolev 
space of order j over R", with i?°(R"^) = L^(R"^). Similarly, we denote by = 
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iy^(R") the weighted Sobolev space defined by the norm = ||a«)||j7j. We also 

set = H^. We write for the intersection fl equipped with the norm 
ll'^llxj ~ II ^11 //J + ll'^^llLj) product X X^~^ equipped with the (e- 

dependent) norm 

\\{wi,W2)\\lj = llwill^,- +e||«;2||x.-i • (1-7) 

Finally, we define = x and Y^a = x L^, equipped with the (e-dependent) 
norms 

||(w^i,w^2)||r, = Ikilllfi +e|k2||i2 , ||(w;i,'u;2)||r,„ = \\iawi,aw2)\\Y, ■ 

Remark that = Y^r\ Y^a- 

It follows from these definitions that {w, wt) G Zl if and only if (atu, awt){'\- + e^^) e 
X L^. Therefore, our perturbation space {{aw^awt) \{w,Wt) e Zl} depends on 7 
and becomes smaller when 7 is increased. On the other hand, using a direct calculation 
in Fourier space, it is not difficult to verify that the origin in (1.6) is linearly stable 
in zl if and only if 1 — 07 + 7^ < 0. In fact, this condition can be read off from the 
coefficient of w in (1.6). As a consequence, the biggest perturbation space in which we 
can hope for stability of the wave is obtained by taking 




Note that this value corresponds to the exponential decay rate of g{x) as a; — > +00, 
since g{x) ~ e"'''^ if c > 2 and g{x) ~ xe~^ if c = 2 [AW]. In the sequel, we shall always 
assume that (1.8) holds, so that (1.6) becomes 

ewu + (1 + 2ec^)wt - 2ecWxt = w^x + ^yW + y/ - 4:Wx - 2gw - auP' . (1.9) 

Furthermore, we shall assume without loss of generality that (7(0) — 1 — a for some 
c < 1/8, and that g{x) > 2a{x)/3 for all x > 0. This can always be achieved by 
replacing g{x) by g{x — xq) for some sufficiently large xq > 0. 

Remark. As in the parabolic case, one can show that the origin in (1.6) is exponentially 
stable in Z^ if c > 2 and 1 — C7 + 7^ < 0. The fastest decay rate is obtained for the 
value 
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Since these results are rather straightforward to prove, we shall focus here on the 
marginal choice (1.8) for which no exponential decay is expected. 

Using these definitions, we can state our first result, which shows that the travelling 
waves are locally stable. 

Theorem 1.1. Assume that n < 4, and let eo > 0, c > 2. Then there exist constants 
Sq > and Kq > 1 such that, for aUO < e < eo, the foUowing holds : for all {^po, ^pi) G 
such that II {(po, 'Pi)\\zi ^ ^^O; there exists a unique solution [w, Wf) G C^(R+, Z^) of (1.9) 
with initial data (to (0), 10^(0)) = ((^05<^i). Moreover, one has 

\\{w{t),wtm\zi < Ko\\{^o,^i)\\zi , (1-11) 

for all t >0, and 

lim {\\Vw{t)\\xo + \\wtmxo + \\w{t)\\L2) =0. (1.12) 
In addition, if c = 2, one has 

lim Vt{\\Vw{t)\\xo + \\wt{t)\\xo + 1^(011^0 = • (1-13) 

t— »+oo ^ " 

Remarks. 

1. ) By a solution of (1.9), we always mean a mild solution, that is a solution of the 
integral equation associated with (1.9), see the proof of Proposition 2.1 below. In gen- 
eral, such solutions satisfy (1.9) in a distributional sense only, see Lions [Li], Section 1.1. 
Remark that Wtt or w^x belong to C°(R+, X~^) only, but Wtt — 2ecWxt — w^x — G 
C°(R+,X°) by (1.9). Moreover, if ((^o^V'i) £ -^e, then the solution {w,wt) belongs to 
C-'^(R+, Zl) n C°(R+, Z^) and satisfies (1.9) in a classical sense. In (1.11) and in the 
sequel, we use the short notation w{t) for ?/;(•,•, t), when no confusion is possible. 

2. ) The restriction n < 4 arises because we control the nonlinearity aw'^ in (1.9) in the 
energy space Z^, using the Sobolev embedding of if^(R") into L'^(R"). More generally, 
if f{u) in (1.1) is a polynomial of degree p > 1, we assume that n < 2p/{p — 1). This 
bound could be improved up to 2{p + l)/(p — 1) using the more sophisticated Lp — Lpi 
estimates of Strichartz [St] , [Br] . 

3. ) If n > 3, it follows from (1.11), (1.12) and the Sobolev embedding theorem that 

2n 

lim {\\w{t)\\L. + \\aw{t)\\L'^) = , 2<q< . (1.14) 

t— *-+oo n — z 
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If n = 2, then (1.14) is valid for all g > 2 and even for = oo if n = 1. In the case 
c = 2, (1.11) and (1.13) imply that 



lim (f^Wwit^L.+t^^^Wawit^L.) =0 

— >+oo \ / 



n{q - 2) 

V = 



4q ' 

for the same values of q. 

Theorem 1.1 is a local result in the sense that the size of the basin of attraction of 
the wave, in particular its dependence on the parameter e > 0, is not specified. However, 
in the parabolic limit e ^ 0, it is known [KR] that the travelling fronts are stable with 
respect to large positive perturbations, and a similar phenomenon is expected to hold for 
(1.5) if e is sufficiently small. To investigate this, we restrict ourselves for convenience 
to one space dimension, and we apply the Maximum Principle for hyperbolic equations, 
which is briefly recalled in Appendix A. Our second result reads: 

Theorem 1.2. Assume that n = 1, and let eo > 0, c > 2, d e (0, 1]. Then for any 
< e < eo and for any constant K > such that 

l-4e(d + K) > , (1.15) 

there exists a constant K* = -fC*(eo, c, d,K)>0 such that the following holds: for any 
{(po,(pi) e Zl satisfying \\{(po,(pi)\\zi < K* and, for (almost) every x e R, 

^o(^) > -{l-d)a{x)-^g{x) , (1.16) 

e(pi{x) > ecipQ{x) - ^{ifio + (1 - d)a~'^g){x) + ec(--fifio + (1 - d)a~'^g'){x) , (1.17) 

there exists a unique solution {w,Wt) G C^(R+, Z^) of (1.9) with initial data {(po, (pi). 
Moreover, one has 

\\{w{t),wtm\zi < K , w{x,t) > -{l-d)a-\x)g{x) , (1.18) 
for allx e'R,te K+, and (1.12), (1.13) hold. 



Remarks. 

1) The proof will show that the constant K* can be chosen so as to satisfy the equation 
KqK*{1 + i^*)i/2 = K/2 for some = KQ{eo,c,d) > 1. Therefore, if e is small, K 
(hence K*) can be chosen very big by (1.15), and Theorem 1.2 shows in this case that 



- 7- 



the travelling wave is stable with respect to large perturbations, provided they satisfy 
the "positivity conditions" (1.16), (1.17). Conversely, if e is large, then K (hence K*) 
has to be very small, and Theorem 1.2 reduces to a local stability result similar to 
Theorem 1.1. 

2) The conditions (1.16), (1.17) appear when applying the Maximum Principle to the 
equation (1.5), see Appendix A. The first one simply says that v(a;,0) > dg{x) for all 
a; e R. The condition on the derivative is not very restrictive if e is small, and disappears 
in the limit e — > 0. 

The previous results are incomplete in the sense that they fail to give a decay rate 
for the perturbations when c > 2. Also, it would be very natural to have at least a 
global existence result if d = 0. Indeed, it is known that, if < v{x,Q) < g{x), the 
solution v{x, t) of the parabolic equation (1.5) with e = exists for all times and satisfies 
< v{x,t) < g{x). In terms of the variable tu, this corresponds to —a~^{x)g{x) < 
w{x,t) < 0. A similar property is expected to hold for the hyperbolic equation (1.9) if 
e is sufficiently small. 

A partial answer to these two questions can be given when e < 1/4. Indeed, in 
this case the Maximum Principle allows us to compare the solution 'w{x,t) of (1.9) with 
solutions of linear equations, whose initial data are the "positive" and "negative" parts 
((/?^, (/?^) of ((/Jo, V'l), ill the sense of Appendix A. They are given by 
^q{x) = sup(0,^o(a;)), 

1 1 (1-19) 

(p"[{x) = c{(p^)'{x) - + ci)(p'^{x) + sup(0, {ipi - cip'q + + ci)^o){x)) , 

and 

V7(7(x) = inf(0, ^o{x)) , 

1 _ 1 (1-20) 

^1 {x) = c(^o y{^) ~ + ^^)^o i^) + inf(0, (^1 - c^o + + c7)^o)(a;)) • 

Remark that ((/?^, ^pf) belong to Z] and that (pi = (p'f -\- (p~ for z = 0, 1. Although the 
norms of (pf seem to depend strongly on e, it is not the case actually: the reader may 
check that |</?^(a;)| < |</?i(a;)| +c|(^Q(a;)| a.e. in R. 

With these definitions, we can state the last result: 

Theorem 1.3. Assume that n = l, and let c>2, d e [0, 1]. Then for any < e < 1/4 
and for any nonnegative constant K satisfying 

1 - 4e(l + K) > , 
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there exists a nonnegative constant K = K{c, K) such that the following holds: for any 
{(fo,(fi) e Zl satisfying (1.16), (1-17) and 

inf(||((^o,<^i)||zi, ||(^^,^^)IUi) <K, 

where {^Pq , (Pi) is given by (1.19), there exists a unique solution (w, Wt) G C°(R+, Z^) 
of (1-9) with initial data {ipo, V'l)- Moreover, one has 

-{l-d)g{x) < a{x)w{x,t) < K , 

for all X e'R, t E R+. Finally, ifd>0 and e < 1/4, one has 

lim t'/^{\\w{t)U^ + \\{w{t),w,mYj = 0. 

t— >+oo 

Remark. The constant K is given hy K = K/N, where N = N{c) is a positive constant. 
Note that the case K = is non trivial: it corresponds to nonpositive initial data. 

An outline of the paper is as follows. In Section 2 we introduce energy functionals 
which allow us to derive a priori estimates for the solutions w{x,t) of (1.9) under the 
assumption that either ||w(t)||xo is sufficiently small or t) satisfies the lower bound 
in (1.18) on some time interval. Using these energy estimates, we prove Theorem 1.1 
in Section 3. Section 4 is devoted to the one-dimensional case n = 1. Combining the 
Maximum Principle with the estimates of Section 2, we derive Theorem 1.2. Further- 
more, when e < 1/4, we obtain linear bounds for the solutions of (1.9) which allow 
us to prove Theorem 1.3. In Section 5, we consider the limiting parabolic equation 
(1.1) when e = 0. Noting that all the estimates made in Section 2 are uniform in e 
when < e < eo, and using the Maximum Principle for parabolic equations, we obtain 
analogues of Theorem 1.1 and Theorem 1.3. Thereby we recover some known stability 
results for the KPP equation. Finally, in Appendix A, we recall the Maximum Principle 
for hyperbolic equations [PW] in a version adapted to our purposes. 



- 9 - 



2. Energy Estimates 

In this section, we derive some a priori estimates for the solutions of (1.9) which will be 
needed in the proofs of Theorem 1.1 and Theorem 1.2. We begin with a standard local 
existence result. 

Proposition 2.1. Let e > 0, c > 2, and let {(po,(pi) e Zl- Then there exists a time 
T = T{e,c,(po,(fii) > such that (1.9) has a unique solution {w,Wt) G C^{[0,T], Zj) 
satisfying {w{0),wt{0)) = {(po,(pi). 

Remark. In fact, the proof gives a lower bound on the existence time which depends 
only on e, c and \\{(po, v'l)!!^!. Moreover, the energy estimates below will show that this 
time is independent of e if e e (0, eo]. 

Proof. Setting W = {w^WtY (where * denotes the transposition), we rewrite (1.9) into 
the "abstract form" 

W = AW + F{W), (2.1) 

where A is the linear operator 

A=( ' ' 

V {dl + Ay + V^^d^ - 2g) -e'^ + 2c(a^ - 7) 

and F{W) = (0, —€~^aw'^Y . It is not difficult to show that the operator A, defined 
on the domain D(A) = Z'^, is the generator of a C*^-semigroup [Pa] of bounded linear 
operators in Z^. Indeed, A can be written as the sum of a bounded operator (depending 
on X E H through the function g) and an unbounded operator with constant coefficients, 
for which the property of being a generator can be verified by a direct calculation 
(using Fourier transforms). Therefore, it follows from a classical stability theorem ([Pa], 
Theorem 3.1.1) that A is the generator of a C°-semigroup e^^ in Z^. 

On the other hand, it is easy to verify that F : — > is a map. Indeed, 
if w G = n H^, then by the Sobolev embedding theorem ||w|||4 < -ftTsHwUl^i 
and ||aw|||4 < i^sHzwH^i for some Ks > 0. Therefore, ||aw^|||2 < Ik'«'|li4||'it'|li4 < 
i^rlll-ujll^i and ||at«^||^2 = ||ow||^4 < Combining these inequalities, 

we find that ||a«)^||x» < -^sll^llxi? which proves that F maps Zj into itself. Since F 
is quadratic, the differentiability follows by the same estimates. 
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In view of these properties, a standard result in semigroup theory ([Pa], Theo- 
rem 6.1.4) imphes that, for aU $ = {ipQ, ipi) e , the integral equation 



JQ 

has a unique solution W e C°([0, T], Z^^) for some T > 0. This is what we call a 
(mild) solution of (2.1), hence of (1.9). Moreover, if $ e Z^, then W e C^([0, T],Zl) n 



In the sequel, we fix eg > 0, c > 2, and for some e G (0, cq] we assume that we are 
given a solution W = {w, wt) of (1.9) (in the sense of Proposition 2.1) defined on some 
time interval [0, T] and satisfying one of the following two assumptions: 

Hypothesis HI: 



Hypothesis H2: 

w{x, y, t) > —(1 — d)a{x)~^g{x) a.e.(x, y),'it e [0, T], for some d G (0, 1]. 

These two cases are adapted to the purposes of the proofs of Theorem 1.1 and 
Theorem 1.2 respectively. To be specific, we assume in the first case that 5 < l/{8Ks), 
where Ks is the constant of the Sobolev embedding of into L'^ (like in the proof of 
Proposition 2.1). 

Under these assumptions, we shall study two families of energy functionals: un- 
weighted and weighted ones, which control the size of the solution w{x,y,t) in the 
spaces and Y^a respectively. We shall derive differential inequalities for these func- 
tionals, which will show that the solution w{x,y,t) is bounded uniformly in time by a 
quantity depending only on the initial data. 

2.1. Unweighted Functionals 

Given w{x, y, t) as above, we define 




C°([0, T], Z^) and satisfies (2.1) in a classical sense ([Pa], Theorem 6.1.5). 



□ 



sup{||t(;(t)||xo 1 1 e [0, T]} < 6 for some sufficiently small 6 > 0, 




(2.2) 



E2{t) 



aEo{t)+Ei{t) , 
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where a = max(2e, l/(2c^)). Here and in the sequel, we set z = {x,y) e and 
dz = dxdy. 



Lemma 2.2. Assume that HI or H2 holds. Then 

Eo{t) > J^^ (^^w^ + + ^9^^^ dz , 

for all t e [0,T]. Moreover, Eq e C^{[0,T]) and 

Eoit) < (c'-4)Eo(t) . 



(2.3) 



(2.4) 



Proof. We first control the cubic term in (2.2). Using the Cauchy-Schwarz and Sobolev 
inequalities, we have 



/ 

Jx<0 



aw^ dz 



< I / a'^w^dz 

'x<0 



1/2 



w'^ dz 



x<0 



1/2 



< 



-^sll'f^^llxo / {u!^ + |V;zt(; p) dz 

Jx<0 



< 2Ks\\w\\x^ 



since g{x) > 1/2 for x < 0. Similarly, 



aw^ dz 



x>0 



< I / w^dz 

'x>Q 



x<0 



1/2 



gw + -I'^zw] ] dz , 



a'^w^ dz 



x>0 



1/2 



< i^5||ti;||xo / I aw 

'x>0 



I (aw^ + \Vz{a^/^w)\A dz . 

Jx>0 ^ ' 



The integral of | V2(a^/^w) |^/^ is equal to 



a;>0 



I i2 7 2 

alV^wj -\- —aw —'^awwx] dz < 



x>0 



(2.5) 



^2 2 ^1 T— 7 I 2 \ 7 

-7 aw H — a VzW dz 
4 ' 2 ' ' ' 



Since 7^ < 1 and 2a{x)/3 < g{x) < 1 for a; > 0, we thus have 



/ 

Jx>0 



aw^ dz 



< 3Ks\\w\\xo 



x>0 



x>0 



< Ks\\w\\xo / ( 2aw'^ H — alV^tyP ) dz 



gw + -\Vzw\ ] dz . 



(2.6) 
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Combining (2.5) and (2.6), we conclude 



aw^ dz 



< 3Ks\\w\\xo 



gw + - \Vzw\ 



dz 



(2.7) 



If HI holds, one has 3K5||w;[|x'' < 3KsS < 1/2, and (2.3) follows immediately. If H2 
holds, then aw^ > —gw^ a.e.(a;,j/) and (2.3) is obvious. 

To prove (2.4), we first assume that {(po,(pi) = {w{0),wt{0)) e Z^. In this case, 
one has {w,wt) G C^{[0,T], Zl), so that Eq e C^([0,T]), and a direct calculation shows 
that 

Eo{t) = -{l + 2ec'y) wldz^^c^ - 4 / w^wtdz . (2.8) 

In the general case where {ipo, (pi) e Z^, we use the fact that the solution (w, wt) G Z^ 
depends continuously on the initial data {(fQ,ipi), uniformly in t G [0,T]. Therefore, 
if G{t) denotes the right-hand side of (2.8), we see that (for fixed t) both EQ{t) and 
Eo{0) + Jl^G{s)ds are continuous functions of {(pQ,(pi) e Z^. Since they coincide on 
a dense subset (namely, Z^), they must be equal everywhere. This proves that Eq e 
C^{[0,T]) and satisfies (2.8). Finally, since 



/ w^Wt dz 






Jr" 





-4 

we see that (2.4) follows from (2.3) and (2.8). 



dz , 



□ 



Lemma 2.3. Assume that HI or H2 holds. Then there exist constants Ki,K2 > 
depending only on eo,c such that 



for all t e [0,T]. Moreover, E2{t) > aEo{t)/2, E2 G C^{[0,T]) and 



(2.9) 



(2.10) 



Remark. We recall that W — {w,wt). The fact that Ki, K2 are independent of e will 
be very important in Section 5, where the limiting case e = is considered. Note that 
the standard choice a = 2e in (2.2) would lead to a constant Ki of order e, see the proof 
below. 
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Proof. Since a > 2e, we have 



\ewwt\ < -w + -€ Wf < -w -\ — aeWf . (2.11) 
' '-3 4^-3 8* ^ ^ 



2 , "ivy 12 , « 2 1 2 1 > 



dz 



Therefore, using (2.3), we find 

E2it) > [ ( 
Furthermore, using (2.7), we obtain 

E2{t) < / (aewt + ^\Vzwf + agw'^ + (1 + ec-f)w'^ + ^aw^) dz 
JR" V 2 3 / 

< (1 + i^5||w||xo) / (aeWf-\ — \V zw\'^ + agw'^ + {1 + ec'y)w'^j 

\ 2 / 

Since < ^ < 1, e < eo, and l/(2c^) < a < max(2eo, l/(2c^)), we arrive at (2.9), with 
Ki,K2 independent of e. Similarly, since \ewwt\ < w"^ /2 + aeWf /A, it follows from (2.2), 
(2.3) that 

Ei{t) > -— j^ytdz > --Eo{t) , 

hence E2{t) > aEo{t)/2. 

To prove (2.10), we proceed along the same lines as in the preceding lemma. Using 
a direct calculation and a density argument, we show that 

Ei{t) = / [ew^ - 2eWxWt - jV^tup - 2gnP - aw^) dz , 

hence 

E2{t) = (e - ci(l + 2ec7)) / wl dz + (a\/ - 4 - 2ec) / w^Wtdz 

Jr" Jr" 

- / (iVztup + 2gw'^ + aw^) dz . 
If Q! = 2e, then ay/ — A — 2ec = — 4e7, and 

I I 2 2 2 2^1 1 2 

Ae'-^\wxWt\ < 4e C7t(;^ H — tu^ < 2eQ!C7t(;^ + -IV^tuI , 



since 7/c < 1/2 by (1.8). If a = 1/(20^) > 2e, then \aVc^ - 4 - 2ec| < ac and 

Cto 22 l|r^i2 

Q!C|Wa;Wt| < -rWf +ac = —Wf +-\yz'W\ . 
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In both cases, we find 

E2{t) ^ - J (^^Wf + zw\'^ + 2gw'^ + aw^^ dz. 

If HI holds, then by (2.2), (2.7) 

1 f f \. 3 5 \ 

E2{t) + -Eo{t) < -J^ {-l\^zwf + -gw^+-aw^] dz 

< -{l-5Ks\\w\\xo) Q|V.«;|' + ^^w;'^ dz < . 

If H2 holds, then we simply have 

E2{t) + Eo{t) < - [ (gw^ + '^awA dz < -\ [ gw^ dz < . 

In both cases, we obtain (2.10). □ 

Remark. Up to now, we did not use the fact that d > 0. Therefore, Lemma 2.2 
and Lemma 2.3 are still valid if H2 holds with d = 0, and the constants Ki,K2 are 
independent of d. 

2.2. Weighted Functionals 

Under the same assumptions as above, we define the weighted functionals 

Foit) = (l^^^t + I'^'^l'^M^ + «W + ^a^w^^ dz , 

F,{t) = J (l^:^aV + ea'^ww^ dz , (2-12) 

F2{t) = aFo{t) + Fi{t) + pEo{t) , 

where a = max(2e, d/{2c^)) and /? = So:. In the case where HI holds, we set d = 1, so 
that a = a. 

Remark. The additional term (3EQ{t) in (2.12) guarantees that F2{t) > 0. However, 
if e is sufficiently small, then aFQ{t) + Fi{t) is already positive, so we may set (3 = 0. 
This possibility will be used in Section 4.3 below. 
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Lemma 2.4. Assume that HI or H2 holds. Then there exist constants K3, K4, K5 > 
such that 

Ks\\Wit)fy^^ < F^it) < K4Wit)fyJl + \\w{t)\\xo)+pEo{t) , (2.13) 
for aU t e [0,T]. Moreover, F2 G C^{[0,T]) and satisGes 

Mt) + ^^F2{t) < K^Eoit) , (2.14) 

where k = d/{8{l + a)). 

Remark. Here and in the sequel, K3, K4, . . . denote positive constants depending only 
on eo, c and, if H2 holds, on d > 0. 

Proof. Using the identity 

/ a^\Vzw\'^dz = I {\V^{aw)\'^ + 'y'^a^w^) dz , (2.15) 

together with the relation 1 — 07 + 7^ = 0, we write F2{t) as 

F2{t) = I ( ya'w^t + ||Vz(a«;)|2 + c7(a/2 - e)a''w^ + |(2^ - l)a''w'' 

a 1 \ ^^-^^^ 

+ —o?w^ + -a?"w^ + ea^wwt ) dz + PEo{t) . 

o A J 

To prove (2.13), we first note that c-^{a/2 — e) < c-^a/2 < a and g <1. Using (2.11), 
we thus find 

F2{t) < j (^aeaV + ||V^(at(;)|^ + (l + a)a2t(;^ + |a^7/^^ d^ + /?£;o(t) . (2.17) 

Furthermore, in analogy with (2.5), we have 

/ a^w^ dz < Ks\\w\\xo / [a'^w'^ + \V z{(iw)\ 2) dz . (2.18) 

Therefore, combining (2.17) and (2.18), we easily obtain the upper bound in (2.13). 
To prove the lower bound, we first use (2.11) and the fact that a > 2e. We find 

F2{t) > + f |V.(att;)|' + |(2(7 - l)a^w^ + + ^aV) dz 

+ PEoit) . 

(2.19) 
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If HI holds, we apply (2.18). Since i^5||t(;||xo < 3/4, we obtain 

F2{t) > j (^|aV + f |V^(aw)|2 + a(^-3/4)aW + ^aW^ dz + fiE^it) . 

If a; < 0, then g{x)-'i/A > l-a-3/4 > 0, since d < 1/8. Ifx > 0, then g{x)-3/4 > -1 
and a{x)'^ < a{x) < 3g{x)/2, so that 

/ a'^w^dz < I / gw^dz < 3Eo{t) , (2.20) 
Jx>o 2 j^>Q 

by (2.3). Therefore, since (3 — 3a, we have 

F2{t) > (^yaV' + ^\^z{aw)f + ^aV^ dz . (2.21) 
If H2 holds, we observe that 

a{g - ll2)c?w'^ + -c?w^ > a{2g/3 - l/2)a^w^ a.e.{x, y) . 

Again, we have 2g{x)/3 - 1/2 > (1 - 4(7)/6 > if a; < 0, and 2g{x)/3 - 1/2 > -1 if 
x >0. Therefore, using (2.19) and proceeding as above, we again arrive at (2.21). This 
proves the lower bound in (2.13). 

To prove (2.14), we proceed along the same lines as in the preceding lemmas. Using 
a direct calculation and a density argument, we first show that 



F2it) = (e — a) / a^Wf dz + c(a — 2e) / a^WxWtdz 

JR" JR" 

- / {\V^{aw)\^ + {2g - l)a'^w'^ + a^w^) dz + PEo{t) . 
Jr" 

2^ 



(2.22) 



If a = d/{2c^) > 2e, then 

I/- -I tin \^ |2 

\c{a - 2e)WxWt\ < ca\WxWt\ < —w^ +-\VzW\ . 
Therefore, using (2.15) and the fact that 7^ < 1, we find 

<- (-a'^w'^ + -\Vz{aw)\'^ + {2g-l-d/2)a'^w'^ + a^wA dz + (3Eo{t) . 
JR" v2 2 / 

(2.23) 

If ci; = 2e, then (2.23) follows immediately from (2.22). 



- 17- 



We now combine (2.17) and (2.23). Using (2.4) and the fact that kck < 1/2, we 
easily find 



(2.24) 



F2{t)+KF2{t) < - J Q|V^(a«;)|2 + (2^-l-d/2-Av(l + d))aW^ dz 
- (1 - Ka/3) / a^w^ dz + (3Eo{t) , 

with /3 = I3{k + - 4). If HI holds, we use (2.18) and obtain 

F2{t) + KF2{t) < - {2g- 3/2 - k{1 + a) - Ks\\w\\xo)a'^w^ dz + /3Eo{t) . 

If x < 0, then 2g{x) - 3/2 - k{1 + a) - Ks\\w\\xo > 1/2 -2a- k{1 + a) - Ks6 > 0, 
by assumptions on cr, S, k. If a; > 0, the same quantity is bounded from below by —2. 
Therefore, using (2.20), we find F2 + KF2 < (6 + /3)Eo, which is (2.14). 

If H2 holds, we infer from (2.24) 
F2{t) + KF2{t) < - {{l + d)g-l- d/2 - k{1 + a))a'^w'^ dz + PEo{t) . 

JR" 

Since g{x) — > 1 as a; — —00, there exists < such that g{x) > 1 —d/8 for all x < xa- 
Therefore, if x < Xd, one has {I + d)g - I - d/2 - k{1 + a) > {3d - d'^) /8 - k{1 + a) > 
by assumptions on d, k. If x > Xd, the same quantity is bounded from below by —2, 
and 

/ a^w^dz < [ gw^dz < ^ e'^'^'''' Eo{t) , 

J X>Xj, ^ J X>Xd ^ 

since 3g{x) > Sg{x — Xd) > 2a{x — Xd) = 2e^^'^a{x) for all x > Xd- Combining these 
inequalities, we find F2 + KF2 < {%e~'^'^^'^ + $)Eq^ which is the desired result. □ 

Corollary 2.5. Assume that HI or H2 holds. Then there exists a constant Kq > 1 
such that 

\\wmz^ < Kewwrniziii + \\wmM'/' , (2.25) 

for allte [0,T]. 



Proof. According to Lemma 2.3, we have 

wmh < ^wmis-^ + ih(o)iixo) , (2.26) 
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for all t G [0, T], since E2 is a decreasing function oft. On the other hand, it follows from 
(2.14) and Lemma 2.3 that F2{t) + nF2{i) < k^E2{t), where = 2K^/a. Integrating 
this inequality, we find 

F2{t) < e-^'F2{0) + K5 [ e-^^'-^^E2{T)dT < F2{0) + —E2{0) , 

Jo « 

hence 

< ^^[K4Wm\h. + {Q + K,/n)K2\\Wm\h){^ + \\H0)\\xo) , (2.27) 

by (2.9), (2.13). Combining (2.26) and (2.27), the result follows. □ 

If H2 holds with d = 0, it is no longer possible to bound W{t) uniformly in time as 
in Corollary 2.5, but the energy estimates above still imply that I|VF(t)||zi cannot grow 
faster than an exponential. This result will be useful in Section 4. 

Corollary 2.6. Assume that H2 holds with d = 0. Then there exist constants p > 
and Kr > 1 such that 

\\Wmzi < K7(l + e''*)||W^(0)|Ui(l + |h(0)|Uo)^/^ , (2.28) 

for all te [0,T]. 

Proof. We recall that Lemma 2.2 and Lemma 2.3 still hold if d = 0, see the remark 
at the end of Section 2.1. Furthermore, if we define F2{t) by (2.12) with a = a = 
max(2e, l/(2c^)) and /3 = 3a, then it is easily verified that (2.13) is still valid. However, 
(2.23) has to be replaced by 

Mt) < -J^ (^|aV + ^|V.(a«;)|2 + ((7-3/2)aV^ dz + PEo{t) • 
Therefore, using (2.4), (2.21), we obtain 

F2it) < 9F2{t) + P{c^ - 4:)Eoit) , 



which replaces (2.14). Integrating this inequality and proceeding as in the proof of 
Corollary 2.5, we obtain (2.28), with p = 9/2. □ 
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3. Local Stability 

In this section, we prove Theorem 1.1 using the energy estimates of Section 2. 
Proof of Theorem 1.1. Let 5o,Ko be defined by the relations 

^6^0(1 + 5o)'/' = S/2 , Ko = Kg{1 + 5o)^/' , 

where Kq > 1 is given in Corollary 2.5 and 5 = 1/(8X5) in the assumption HI, 
Section 2. Then, for all ^> = {(fo,(fi) e Z\ such that ||^||zi < (^o, Eq.(1.9) has 
a unique global solution Wit) = {w{t),Wt{t)) e satisfying W{0) = Indeed, 
in view of the local existence result (Proposition 2.1), it is sufficient to show that 
||W^(^)IUi < <^ whenever W{t) exists. Assume on the contrary that there exists a 
time T > such that ||W^(T)||^i = S and ||W^(t)|Ui < S for all t e [0,T). Then 
lk(t)||xo < ||W^(*)IUi < S for all t e [0,T], so that HI holds on [0,T]. By Corol- 
lary 2.5, it follows that ||W^(T)||^i < KeSoil + 5o)^/^ = 5/2, which is a contradiction. 
Therefore, W{t) exists for all times and the assumption HI is always satisfied. By 
Corollary 2.5 again, we conclude that ||VF(t)||zi < -f^oH^IUi for all t > 0, which proves 
(1.11). 

To prove (1.12), (1.13), we use the differential inequalities satisfied by the func- 
tionals i?o,-E'2,-^2 defined in Section 2. The following arguments are standard (see for 
example [EW]) and will be reproduced here for the sake of completeness. First, since E2 
is a positive, decreasing function of t by Lemma 2.3, E2{t) converges to a nonnegative 
limit as t +00. By (2.4), (2.10), so does Eo + 2{c^ -A)E2. Therefore, Eo{t) converges 
as t — > +00, and since i?o(^) > it follows from (2.10) that the integral 



I 



+00 

Eo{T)dT < 2(^2(0) -E2(+00)) 



is finite. Thus Eo{t) converges to zero as t — > +00. Moreover, integrating the differential 
inequality (2.14), we find 

F2it) < e-''*F2(0) + K5 f e-^^'-^^EQ{T)dT 

Jo 

< e-''*F2(0) + K5 (e-'^*/' f^'' Eo{T)dT+ f e'^^'-^^ Eo{r) dr] (3.1) 

\ ^0 Jt/2 J 

< e-''*/2(F2(0) + 2X5^^2(0)) + — sup Eo{t) , 
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hence F2{t) converges to zero as t — > +00. Therefore, using the lower bounds in (2.3), 
(2.13), we obtain (1.12). 

In the case where c = 2, Eq itself is a decreasing function of t by (2.4), hence 
tEo(t) < 2 J^^^Eo(T)dT. Thus tEo(t) converges to zero as t — > +00, and by (3.1) the 
same is true for tF2{t). Therefore, using again (2.3), (2.13), we obtain (1.13). This 



4. Global Stability Results in the One-Dimensional Case 

Throughout this section, we assume that n — 1. First, we prove Theorem 1.2 using 
the results of Section 2 and the Maximum Principle for hyperbolic equations. Then, we 
study in more details the case e < 1/4; we give linear upper and lower bounds for the 
solutions of (1.9). Finally, using these linear bounds, we prove Theorem 1.3. 

4.1. Global Stability in the General Case 

We first show that the assumption H2 (see Section 2) holds if the solution w{x,t) of 
(1.9) is bounded from above and if the initial data satisfy (1.16), (1-17). 

Proposition 4.1. Let e > 0, (i e [0, 1] and let K be a nonnegative constant such that 



For some T > 0, assume that {w,Wt) G C°([0,T], Zl) is a solution of (1.9) with initial 
data {(po,(pi) satisfying (1.16), (1.17), namely 



concludes the proof of Theorem 1.1. 



□ 



1 - 4e(d + K)>0 . 



(4.1) 



Vo{x) > -{l-d)a ^{x)g{x) , 



(4.2) 



e(pi{x) > ec(fiQ{x) - -{ipo + {l- d)a ^g){x) + ec{--f(po + (1 - d)a ^g'){x) , (4.3) 



for (almost) every x eH. Suppose moreover that 



a{x)w{x,t) < K , V(a;,t) e R X [0,T] . 



(4.4) 



Then 



w{x,t) > -(1 - d)a-^{x)g{x) , V(a;,t) e R x [0, T] . 



(4.5) 
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Proof. We recall that the inequality (4.5) is equivalent to v{x,t) > dg{x), where 
v{x,t) is the solution of (1.5). Also, we remark that dg — v belongs to the space 
C0([0, T],Hl^{-R)) n CH[0, T],Ll^{-R)) and satisfies 

L{dg — v) = C{dg — v) + h{x, t){dg — v) = g^d{l — d) , 

where £(f) = v^x + '^^cvxt — ^vu + cvj^ — vt a,ndh{x,t) = l — v{x,t)—dg{x). Therefore, 
to prove (4.5), we are led to apply the Maximum Principle (Theorem A.l, Appendix A) 
to the function dg — v and to the operator L. The condition (A.l) is obviously satisfied. 
Due to (4.4), the condition (A. 2) holds, i.e., 

h{x, t) > 1 - (1 + d)g{x) -K>-K-d. 

This estimate and (4.1) imply that (A. 3) also holds. Moreover, since < (i < 1, we have 
L{dg~-v){x,t) > Oa,.e.{x,t) G Rx[0,T], which is (A. 4). Finally, the conditions (A. 5) and 
(A. 6) required on dg — v are nothing else but the hypotheses (4.2) and (4.3). Therefore, it 
follows from Theorem A.l that {dg — v){x,t) < 0, that is, w{x,t) > —{1 — d)a{x)~^ g{x) , 
for aU (a;,t) e R X [0,T]. □ 

Using Proposition 4.1 and Corollary 2.5, we now prove the first global stability 
result. 

Proof of Theorem 1.2. The proof is very similar to the one of Theorem 1.1 given in 
Section 3. Let be a real number, < < 1. We define K* by the relation 

where > 1 has been introduced in Corollary 2.5. According to Proposition 2.1, 
there exist a time T > and a unique solution W{t) = {w, Wt) G C°([0, T), Zj) of (1.9) 
with initial data ((/7o,9?i) such that ||M^(t)||zi < K for all t G [0, T) and, if T < cxd, 
W{t) G CO([0,T], Zl) and ||W^(T)||^i = K. We show by contradiction that T = oo. If 
T < cxD, we have 

a{x)w{x,t) < \\aw(t)\\Loo < \\aw(t)\\Hi < \\W{t)\\zi < K , 

for all t G [0,T]. Thus, by Proposition 4.1, w{x,t) > —(1 — d)a{x)~^g{x) for all 
{x,t) G R X [0, T], i.e. the assumption H2 of Section 2 holds on [0, T] (we recall that 
d > Ohere). By Corollary 2.5, it follows that \\W{T)\\zi < Kq{1+K*Y/^K* = {l-^i)K, 
which is a contradiction. Therefore T = oo, and the inequalities (1-18) hold for all times. 
The properties (1-12), (1.13) are proved like in Section 3. □ 
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4.2. Linear Bounds in the Case e < 1/4 

From now on, we assume that e < 1/4. In this case, the range of apphcation of the 
Maximum Principle is much wider, and we can show that the solution w{x,t) of (1.9) 
is bounded from above and from below by solutions of suitable linear equations. These 
linear bounds will be crucial for the proof of Theorem 1.3. Before stating the results, 
we introduce some additional notation. 

For all d G [—1, 1], we denote by S'd(t) e ^^{2^, Z^) the linear group associated with 
the equation {L(iw){x,t) = 0, where 

LdW = -ewtt - (1 + 2ec^)wt + 2ecwxt + Wxx + a/c^ - 4 tUa; - (1 + d)gw . (4.6) 
For {(fQ, ipi) e , we set 

-Sd(t)(^o,^i) = {wd{t),Wdt{t)) . (4.7) 

In (1.19), (1.20), we have defined the positive and negative parts {^^^^^i) of ((/7o,</?i)- 
In analogy with (4.7), we set 

Sdmvtvt) = {w^{t),w%{t)) . (4.8) 

We now show the existence of a linear upper bound. 

Lemma 4.2. Let e < 1/4. For any {(fQ,<pi) G Zl, the solution {w,Wt) e C^{[0,T],Zl) 
of (1-9), with initial data {(fio,(pi), satisGes for any d E [—1,1], 

w{x,t) < wi{x,t) < wj{x,t) , y{x,t)eRx [0,r] , (4.9) 

where wi, w'^ have been defined in (4.7) and (4.8) respectively. 

Proof. We first prove the inequality w{x,t) < wi{x,t). The function w — wi satisfies 
the equation Li(w — wi) = aiu^ > 0, where Li has been defined in (4.6). Thus, we can 
apply the Maximum Principle to the function w — wi and to the operator Li. Indeed 
the conditions (A.l), (A. 2), (A. 4) are satisfied, and, since the initial data for w,w 
coincide, (A. 5) and (A. 6) obviously hold. Since —2g{x) > —2, the condition (A. 3) with 
h = —2 becomes (1 — 4e)(e + e^c^) > 0, which is satisfied because e < 1/4. Therefore, 
Theorem A.l. implies that w{x,t) - wi{x,t) < for all {x,t) e R x [0,T]. 



- 23 - 



We next show that w'^{x,t) > 0. Since Ld{—w'^) = 0, we can apply the Maximum 
Principle to the function —w'^ and to the operator L^- In view of the first part of 
the proof, the conditions (A.l) to (A. 4) hold. Due to the choice of ((/j^J", t/?^) made in 
(1.19), the hypotheses (A. 5) and (A. 6) are also satisfied. Therefore, by Theorem A.l, 
w'^ix, t)>0 for aU {x, t) e R x R+. 

Finally, we show that wi{x,t) < w^{x,t) for all d e [—1,1], by applying the 
Maximum Principle to the function wi — w'^ and to the operator Li. As we have 
already remarked, the hypotheses (A.l), (A. 2), (A. 3) are satisfied. The condition (A. 4) 
holds, since Li{wi — w^) = (1 — d)gw'^ and > 0. The choice of {(p^, (pf) in (1.19) 
also implies that (A. 5) and (A. 6) hold. Hence, by Theorem A.l, wi{x,t) < w'^{x,t) for 
all (a;,t) e R X [0,T]. □ 

In a similar way, we obtain linear lower bounds for w{x,t). 

Lemma 4.3. Let e < l/A, d E [0, 1], and let K be a nonnegative constant such that 

1 -4e(l + K) > . (4.10) 

For sonic T > 0, assume that (to, Wt) G C^{[0, T], Z^) is a solution of (1-9) with initial 
data {(fio, (fii) e satisfying (4.2) and (4.3). Suppose moreover that (4.4) holds. Then 

wZi{x,t) <w^(x,i) <w{x,t) , V(a;,t) e R X [0,T] , (4.11) 

where wZi, w'^ have been defined in (4.8). 

Proof. As in the proof of Lemma 4.2, we show that w'^ {x, t) < by applying the Max- 
imum Principle to the function w'^ and to the operator L^. To show that wZi{x,t) < 
{x, t), we apply the Maximum Principle to the function wZi — w'^ and to the operator 
L-i. Since L-i{wZi — w^ ) = —(1 + d)gw~^ and to^ < 0, the hypothesis (A. 4) holds. 
The other conditions are obvious or have been verified in the proof of Lemma 4.2. 

It remains to prove that w(x.,t) > w'^{x,t). We again apply the Maximum Prin- 
ciple, but now to the function w'^ — w and to the operator L| = L_i + h*, where 
h*{x,t) = -{2g{x) + a{x)w{x,t)). Since h*{x,t) > -2 - K, the condition (A.3) be- 
comes (1 — 4e(l -|- K)){e + e^c^) > 0, which is nothing but (4.10). Moreover, we have 
Ll{w~i^ — w) = —w~i^ {aw+{l—d)g) > 0, since w'^{x, t) < and a{x)w{x, t) > —{l—d)g{x) 
by Proposition 4.1. Thus (A. 4) holds, and due to the choice of {(pQ,(fi^) in (1.20) the 
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conditions (A. 5) and (A. 6) are also satisfied. Therefore {x, t) < w{x, t) by Theo- 
rem A.l, and Lemma 4.3 is proved. □ 

Since wi{x,t) is a solution of the linear equation Liw = 0, it is easy to bound it in 
terms of the initial data {(po, (pi). We have the following result: 

Lemma 4.4. Let € < 1/4:. There exists a constant N = N{c) > 1 such that 

||'5i(t)|U(zi,zi) <A^, VteR+. (4.12) 

Proof. All we need is repeating the energy estimates of Section 2 for the linear equation 
obtained by dropping the last term —aw'^ in (1.9). The functionals Eo{t) , E2{t) , F2{t) 
are then replaced by the quadratic expressions 

Eo{t) = + ^h^r + ^«^') dx , E2{t) = aEo{t) + E^{t) , 

hit) ^ a j [^a^wl + ^aVxl^ + dx + F^it) + /3^o(t) , 

where a, /3, £'i(t), Fi(t) are defined in (2.2), (2.12). Of course, the assumptions HI, H2 
are not needed anymore, since they were used to control the cubic terms in E2{t), F2{t). 
Following exactly the lines of the proof of Lemma 2.2, Lemma 2.3, Lemma 2.4 (with obvi- 
ous simplifications), we arrive at Corollary 2.5, which reduces in this case to || VF(t) < 
iV||l¥(0)||zi for some N{c) > 1. This proves (4.12). □ 



4.3. Global Stability and Decay in the Case e < 1/4 

Using the linear bounds of the previous paragraph, we are now able to improve the 
global stability results. Theorem 1.3 will be a direct consequence of the following two 
propositions: 

Proposition 4.5. Let e < l/A, d E [0, 1] and K be a nonnegative constant, such that 
(4.10) holds. For any {(po, cpi) e Z\ satisfying (4.2), (4.3) and 

inf(||(^o,^i)IUi,||(^o+,^^)IUi) < ^ , (4.13) 

where {(fQ,(fi) is defined in (1.19) and N in Lemma 4.4, there exists a unique global 
solution {w,Wt) e C°(R_|_, Z]) of (1.9) with initial data {(po, cpi). Moreover, we have 

-(1 - d)g{x) < a{x)w{x, t) < K , V (x, t) e R x R+ , (4.14) 
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and 

wZi{x,t) < w{x,t) < wi{x,t) < wt^{x,t) , V(x,t)eRxR+ . (4.15) 
In addition, if d > 0, the properties (1-12) and (1-13) hold. 

Remark. The case K = is non trivial, because it corresponds to {<PQ,<pt) = 0, i.e. 

(fio{x) < 0, (fii{x) < c(fQ{x)-{^ + c'y)(fo{x). 

In this case, (4.14) shows that w{x,t) stays nonpositive for aU times. 

Proof. According to Proposition 2.1, there exist a maximal time T > and a solution 
{w,wt) e C°([0, T), Zg^) of (1.9) with initial data {ipQ,(^i) such that either T = cx), or 
T < oo. In the latter case, there exists a sequence of positive times tn, tn < T, such 
that tn ^ T as n ^ +oo and 

\\iw{tn),Wtitn))\\zl ^ +00 , (4.16) 

as n — > +00. By Lemma 4.2 and Lemma 4.4, we have 

a{x)w{x,t) < inl{a{x)wi{x,t), a{x)wi{Xjt)) < inf(||t(;i(t)||xi, ||M'i^(t)||xi) 
< Nmm^o.^i)\\ziAM.V^t)\\zi) < K , 

for all {x,t) e R X [0,T). Therefore, by Proposition 4.1, 

w{x, t) > -(1 - d)a{x)-^g{x) , V (a;, t) e R x [0, T) , 

which allows us to apply Corollary 2.6. It follows that 

\\{w{t),wtm\zi < K7(l + e^^)||(¥'o,¥'i)||zi(l + lbo|Uo)'/' , Vt e [0,T) , (4.17) 

which contradicts (4.16). Thus T = oo. The estimates (4.14), (4.15) are direct conse- 
quences of Proposition 4.1, Lemma 4.2 and Lemma 4.3. If d > 0, the properties (1.12), 
(1.13) are obtained like in the proof of Theorem 1.1 in Section 3. □ 

Remark. If d = 0, we can still show, by arguing as in the proof of Theorem 1.1 in 
Section 3, that 

lim {\\W:,{t)\\L2 + \\wt{t)\\L2) = , 
t— >+oo 
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and that this quantity is 0{t ^/^) if c = 2. However, since (2.14) no longer holds, we 
cannot show that || (^w, '*^t)|| y^a converges to zero in this case. 

Proposition 4.6. Under the assumptions of Proposition 4.5, the solution {w,Wt) G 
C^(R_i_,Z^^) of (1-9) with initial data {ipQ,ipi) satisGes 

lim t^/^||i(;(t)||Loo = . (4.18) 

If, in addition, d> and e < 1/4, then 

hm tV4(||^(^)||^^ + ||(^(^),^^(^))||^^j = 0. (4.19) 

t— »+oo 

Proof. We first prove (4.18). From (4.15), it follows that 

lk(t)||Loo < sup(||t(;li(t)||Loo,||t(;Ii(t)||ioo) , VteR+. (4.20) 

Therefore, we need only show that (4.18) holds for any solution w G C*^(R-|-,X^) fl 
C^(R-l-, X^) of the linear equation with constant coefficients L-iw = 0. Again, this can 
be done using the energy estimates of Section 2. Indeed, setting w{x, t) — u>{x + ut, t), 
where u — \/c^ — 4/ (1 + 2ecj), we see that u{x,t) satisfies 

eutt + (1 + 2ec7)a;t - 250;^;^ = Aco^^ , (4-21) 

where B > and A = (1 + ec^ + ecVc^ - 4) /(I + 2ecj) > 0. Remark that the coefficient 
of (jJx vanishes in (4.21), like for the equation (1.9) in the case c = 2. Therefore, setting 

Eo{t) = {^ojI + ^A|a;^|2^ dx , 

and modifying accordingly the definitions of Ei{t) and E2{t), we show like in Section 2 
that -E'o(^) ^ and that (2.9), (2.10) hold. Then arguing like in Section 3, we find 
that ||a;(t) < Co||a;(0)||//i, and limt^_|_oo ^"'"^^||'^x(^)||l2 = 0. Since ||c<;(t) H^oo < 
ll^^l^) IIl2^II'^x(^)||^2^) we see that lim^_>_|.oo t'''/^||a;(t)||/,oo = 0, which together with (4.20) 
proves (4.18). 

To prove (4.19), we recall that, if e < 1/4, we can define F2{t) by (2.12) with /? = 0. 
Indeed, if e = 1/4 — (5 for some 5 > 0, it is easy to verify that, under the assumption H2, 



F2{t) > SKs{c)\\W{t)\\l^ , 



(4.22) 
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for some constant Ks{c) > 0. Proceeding again like in the proof of Lemma 2.4, we show 
that 

F2{t)+KF2{t) < Ci j a^{x)w^{x,t)dx < C2\\w{t)\\l^ , (4.23) 

Jx>Xd 

where C2 — {Ci/2'y)e~'^'^^'' . Integrating this differential inequality and using (4.18) 
and (4.22), we obtain (4.19). The proof of Proposition 4.6, hence of Theorem 1.3, is 
complete. □ 

Remark. Since L_i is a linear operator with constant coefficients, it is possible to 
obtain explicit expressions for the solutions of the equation L^iw = in terms of the 
initial data, see for example [Sm], chap. VII. 2. 5. Therefore, (4.18) could also be proved 
by a direct (but cumbersome) calculation. 



5. The Limiting Case e = 
5.1. Local Stability 

If we set e = in (1.1), we obtain the well-known parabolic KPP equation, the travelling 
wave solutions g{x) of which are given by (1.2) for c > 2. To study their stability, we 
proceed like in the Introduction. First, using the change of variables (1.4), we arrive at 
(1.5) with e = 0. Then, we look for solutions of the form u(x, y, t) — g{x)+a{x)w{x, y, t), 
where a{x) ~ e~^^ , and we are led to study the stability of the solution = of the 
parabolic equation (1.6) for e = in the Sobolev space = fl H^. Again, linear 
stability holds if and only if 1 — 07+7^ < 0, so the biggest perturbation space is obtained 
by choosing 7 as in (1.8). Then, the equation (1.6) for e = becomes 

wt = Wxx + ^yW + s/ — Awx — 2gw — aw^ . (5-1) 

It is known in this case that the origin is stable in X^, with polynomial decay of the 
perturbations to zero as t — > -|-oo. 

Remark. In the case c > 2, it is also known that the origin is exponentially stable in 
X^ if 1 — C7 + 7^ < 0. The best decay rate is obtained for the value 7 = c/2 [Sa], which 
is precisely (1.10) for e = 0. 

In Section 2, we have introduced various energy functionals for e > 0, which were 
used to estimate the different norms of the solution {w,wt) of (1.9). These functionals 
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are all well defined for e = and allow us to control the norm of the solution w of (5.1). 
Since all the estimates are uniform in e as e goes to 0, we can follow exactly the lines of 
the proof of Theorem 1.1, and we arrive at the (already known) local stability result: 

Theorem 5.1. Assume that n < 4, and c > 2. Then there exist constants 6o > and 
Kq > such that the foUowing holds: for all (fo ^ satisfying \\fo\\x'^ ^ ^Oi there 
exists a unique solution w G C^{H-^-, X'^) of (5.1) with initial condition w{0) — (po. 
Moreover, one has \\w{t)\\x^ < -f^o||</'o||xi for all t > 0, and 

lim (||V«;(t)|Uo + +Mt)|U^) = 0. 

In addition, if c= 2, one has 

hm Vt{\\Vw{t)\\xo + \\w{t)h.) =0. 

Remark. Contrary to the hyperbolic decay rate in time of the solution w{t) 

of (5.1) is easily obtained for all c > 2. Indeed, following the ideas of Nash, it is a 
classical task to estimate the L^-norm of solutions to parabolic equations for p > 2. In 
our case, if we know an upper bound on ||t(;(t)||LP, then we can show that ||t(;(t)||^2p 
decays to zero (like an inverse power of t) as t — > oo, see [FS]. Thus, using the L^-bound 
of Theorem 5.1 and proceeding by recursion, we can show that 

\\w{t)\\L. + \\w{t)\\L, = O (t-^^) , t^+OO, 

where rj — n{q — 2)/(4g) and g > 2 is as in (1.14). 
5.2. Global Stability 

Like in the hyperbolic case, we obtain a global stability result when n = 1. But here we 
apply the Maximum Principle for parabolic equations on unbounded domains as given 
in [PW], Section 3.6. Remark that it is required that w{x,t) does not grow faster than 
exp(Ca;^) as x goes to ±oo (uniformly in t), a condition which is clearly satisfied in 
our case. Like in Paragraph 4.2, we denote by Sd(t) e C{X^,X^) the linear semigroup 
associated with the equation 

wt = + V c2 -Awx - (1 + d)gw , d e [-1, 1] . 
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For ipo G X^, we set 'Ed{t)(po = Wd{t), Ed(t)(/?J = w^(t), where (p^ have been given 
in (1.19), (1.20). Then, following the lines of the proofs contained in Section 4, we 
obtain the global stability result below, which has already been known, though maybe 
not exactly in this form. 

Theorem 5.2. Assume that n = 1, and let c > 2, d E (0, 1]. Then, for any (fo G 
satisfying (1.16), namely ^po{x) > —{1 — d)a{x)~^g{x) for all x G R, there exists a 
unique solution w G C°(R_|_,X^) of (5.1) with initial condition w{0) — ipo- Moreover, 
one has w{x, t) > —(1 — d)a~^{x)g{x), and 

wZi{x,t) < w{x,t) < wi{x,t) < w^i{x,t) , 

for all x e'R, t E R+. In particular, if d > 0, 

lim t'/m\w{t)\\L^ + \\wmHi) =0. 



Appendix A. Maximum Principle for a Hyperbolic Operator 

We consider the following hyperbolic operator C with constant real coefficients 

C{u) = Auxx + 2Buxt + Cutt + Dux + Eut , 

where 

C < , - AC > . (A.l) 
We introduce a function h G C^(R x [0,T]) satisfying 

h{x,t) > h, for aU {x,t) G R x [0, T] , (A.2) 

where T is a positive number and his a real number. We suppose in addition, that the 
condition 

{E'^ - 4Ch){B^ - AC) > {BE - CDf , (A.3) 

holds. Finally, we set L = £ + h{x,t). The following Maximum Principle is a sim- 
ple consequence of the one given by Protter and Weinberger (see [PW], Chapter 4, 
Theorem 1). 
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Theorem A.l. Assume that the conditions (A.l), (A. 2) and (A. 3) hold. If the function 
{u{x, t),ut{x, t)) belongs to C^{[0, T],Hl^{'R) x Ll^{Il)), with Au^a> + '^Bu^t + Cuu in 
L^^^(R X (0, T)), and satisfies the following properties, 

L{u){x, t)>0 , a.e.{x, t) eKx [0, T] , (A.4) 

u(x,0)<0, Va;eR, (A.5) 
-Cut{x,0) - Bua;{x,0) -^Eu{x,0) <0 , a.e. , (A.6) 
then u{x,t) < for all {x,t) e R x [0,T]. 

Proof. Plotter and Weinberger proved their Maximum Principle under the stronger 
assumption u{x,t) e C^(R x (0,T)) fl C^(R x [0,T)), but their proof generahzes 
easily to functions u satisfying the weaker regularity hypothesis {u{x,t),ut{x,t)) G 
C°([0, T],Hl^{n) X LL(R)), with Au^^ + 2Bu^t + Cuu in Ll^{Il x (0, T)). Indeed 
their key identity (see [PW], Equation (3), page 202) still holds under these weaker 
regularity assumptions and is proved by a density argument. 

If E = D = 0, the result of Theorem A.l is a direct consequence of the above 
remark and of Theorem 1 of [PW]. Indeed, thanks to our assumptions (A.l), (A. 2), 
(A. 3), the condition of [PW] on the operator L, that is h{x,t) > 0, is clearly satisfied. 
Since = 0, the conditions required on w(a;, t) are exactly the hypotheses (A.4) to 
(A.6). 

If £^ 7^ or D 7^ 0, we reduce our problem to the previous case by introducing the 
function 

v{x,t) = e-'^^-l^''u{x,t) , 

where 

EA-BD ^ CD- EB 
a = — — - , p — 



2(^2 - AC) ' 2(^2 - AC) ■ 

A short computation shows that 

L{u){x,t) = e«*+^^L(^;)(x,t) , 

where L{v) — Avxx + '2Bvxt + Cvu + hv and 

=-^{(E2_4CMx,t))-^^^^}. 
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Now, we can apply the previous Maximum Principle, where D = E = 0, to the operator 
L and to the function v. Indeed, due to the hypotheses (A.l), (A. 2), (A. 3), h{x,t) > 
for aU {x,t) e R X [0,T]. Moreover, L{v){x,t) > a..e.{x,t) e R x [0, T] and v satisfies 

v{x, 0) < , -Cvt{x, 0) - Bva:{x, 0) < , 

which are exactly the required conditions. Thus, we have proved that v{x, t) < 0, hence 
u{x, t) <0 for aU {x, t) e R x [0, T] . □ 

Remark. Theorem A.l suggests the following definition of a partial order in i/^^^^(R) x 
L^^^(R). We say that {^0,^1) < {i^o,i^i) if 

<fioix) < ipoix) , Va; e R , 
-C(pi{x) - B(p'q{x) - ^E(fio{x) < -CiIji{x) - BiIjq{x) - ^EiIjq{x) a.e. , 

see (A. 5), (A. 6). Then, if {(fo,(fi) < {t(jo,t(ji), the solution of the linear hyperbolic 
equation L{u){x, t) = satisfying u{x, 0) = v?o(a^), Utix, 0) = (pi{x) stays for all t G R+ 
below the corresponding solution with initial data (i/'ojV'i)- An important property of 
this order is that we can write any {<fo,(f>i) G H^^^ x Lf^^ as the sum of a "positive" 
part (ipQ^ipf) > and a "negative" part ((pQ^ipi) < 0. This decomposition is unique 
if we impose that {(pQ^ip^) = whenever (ipo^ipi) < and {ipQ,(p^) = whenever 
{^0, fi) > 0. In the case of the operator Ld defined in (4.6), for which C = — e, B = ec, 
E = —(1 + 2ec7), the corresponding formulae for {(pg, (pf) are given in (1.19), (1.20). 
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